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ABSTRACT 

Observations of H2O masers from circumnuclear disks in active galaxies for the Megamaser Cos- 
mology Project allow accurate measurement of the mass of supermassive black holes (BH) in these 
galaxies. We present the Very Long Baseline Intcrfcromctry (VLBI) images and kinematics of water 
maser emission in six active galaxies: NGC 1194, NGC 2273, NGC 2960 (Mrk 1419), NGC 4388, 
NGC 6264 and NGC 6323. We use the Keplerian rotation curves of these six megamaser galaxies, 
plus a seventh previously published, to determine accurate enclosed masses within the central ~ 0.3 
pc of these galaxies, smaller than the radius of the sphere of influence of the central mass in all cases. 
We also set lower limits to the central mass densities of between 0.12 and 60 x 10 10 M© pc~ 3 . For six of 
the seven disks, the high central densities rule out clusters of stars or stellar remnants as the central 
objects, and this result further supports our assumption that the enclosed mass can be attributed 
predominantly to a supermassive black hole. The seven BHs have masses ranging between 0.76 and 
6.5xl0 7 M Q . The BH mass errors are « 11%, dominated by the uncertainty of the Hubble constant. 
We compare the megamaser BH mass determination with other BH mass measurement techniques. 
The BH mass based on virial estimation in four galaxies is consistent with the megamaser BH mass 
given the latest empirical value of (/), but the virial mass uncertainty is much greater. Circumnuclear 
megamaser disks allow the best mass determination of the central BH mass in external galaxies and 
significantly improve the observational basis at the low-mass end of the M — a* relation. The M — a* 
relation may not be a single, low-scatter power law as originally proposed. MCP observations continue 
and we expect to obtain more maser BH masses in the future. 

Subject headings: accretion, accretion disks - galaxies: nuclei - galaxies: masers - galaxies: active - 
galaxies: ISM - galaxies: Seyfert 



1. INTRODUCTION 

The primary goal of the Megamaser Cosmology Project 
(MCP; Braatz et al. 2009; Reid et al. 2009a; Braatz et al. 
2010) is to determine the Hubble constant Hq to ~ 3% 
accuracy in order to constrain the equation of state pa- 
rameter w of dark energy. The key to achieving this goal 
is to measure accurate distances to galaxies well into the 
Hubble flow (50 — 200 Mpc). A proven method to mea- 
sure accurate angular-diameter distances involves sub- 
milliarcsecond resolution imaging of H2O maser emis- 
sion from sub-parsec circumnuclear disks at the center of 
active galaxies, a tech nique established by t he study of 
NGC 4258 with VLBI ( jHerrnstein et aTJ l999h. Through- 
out the paper, we will refer to such circumnuclear disks as 
"megamaser disks" and galaxies that contain such mega- 
maser disks as "megamaser galaxies". As the MCP dis- 
covers and images megamaser disks, an important result, 
in addition to the distance determination, is the accurate 
measurement of the enclosed masses of the massive dark 
objects (MDOs; Kormendy & Richstone 1995; Richstone 
et al. 1998) at the centers of these megamaser galaxies. 
We assume that the MDOs are black holes (BHs) and 
justify this assumption in section 4. 



The VLBI plays a crucial role in measuring BH masses 
with high precision especially in galaxies with lighter su- 
permassive black holes (i.e. M. ~ 10 6 - 10 7 M Q ). The 
critical advantage provided by VLBI is angular resolution 
two orders of magnitude higher than the best optical res- 
olution. For any given galaxy with a nearly constant cen- 
tral mass density of stars, Mbh is proportional to Rf nf , 
where i?j n f is th e radius of the gr avitational sphere of 
influence of BH (jBarth et al.ll2003f K So, a factor of 100 
increase in resolution permits measurements of masses 
up to 10 6 times smaller. Similarly, the central density 
limits that can be set are up to 10 6 times higher, high 
enough to rule out extremely dense star clusters as the 
MDOs based on dynamical argument (see section 4) in 
most megamaser disks presented here. 

Megamaser disks, such as the archetypal one in NGC 
4258, are usually found in the centers of Seyfert 2 galax- 
ies. Because maser emission is beamed and long path 
lengths are required for strong maser amplification, the 
megamaser disks are observable only if the the disk is 
close to edge-on. In NGC 4258 the disk is inclined —82° 
and its rotation curve is Keplerian to better than 1% 
accuracy, which makes the BH mass determination ro- 
bust with very few assumptions. The megamaser disks 
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are typically smaller (r ^0.2 pc in NGC 4258) than the 
gravitational sphere of influence of their supermassive 
black holes (r ~ 1 pc in NGC 4258); this guarantees that 
the gravitational potential is dominated by the central 
mass. Applying dynamical arguments to the extremely 
high mass density po ~ f 11 M Q pc~ 3 within the NGC 
4258 maser disk can rule out a dense cluster of st ars or 
stella r remnants for most of the central dark mass (|Maozl 
119951) . implying that the measured central mass is dom- 
inated by a BH. 

The megamaser disk method for estimating BH masses 
has some practical limitations. First, finding mega- 
maser disks is difficult, partly because detectable disks 
need to be nearly edge-on. Only eight BH masses have 
been publis hed based on measur ements of megamasers : 
NGC 1068 (|Greenhill et al.l fl996). NGC 2960 (Mrk 1419) 
(Henkel et al. 2002; n ot based on V LBI), NGC 3079 
(IKondratko et al.l [2005] ), NGC 3393 (jKondratko et al.l 
120081), UGC 3789 (Reid et al. 2009 ; MCP pa- 
per I), NGC 4258 dHeTrnstein et al.l [T999h . NGC 4945 
(iGreenhill et al.l 11997ft . and Circinus (jGreenhill et al.l 
120031) . Second, some rotation curves are significantly flat- 
ter th an Keplerian, e.g. NG C 1068; NGC 3079; and IC 
1481 (|Mamvoda et al.1 12009). The origin of the flatter 
rotation curves in these galaxies is unclear. It could be 
caused by self-gravity of a massive disk the presence of a 
nucle ar cluster, or radiation pressure ([Lodato fc Berth] 
2003). In these cases, the measured enclosed mass of 
the disk could be significantly different from the actual 
BH mass. Therefore, the key to measuring reliable BH 
masses using megamasers is to find megamaser disks that 
have well defined edge-on disks with Keplerian rotation 
curves. 

Over the past two decades, there has been substantial 
progress in detecting BHs and constraining their masses, 
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rapid progress was mainly facilitated by the high angu- 
lar resolution provided by the Hubble Space Telescope 
(HST) and by gradually maturing techniques for mod- 
eling stellar dynamics of galaxies. The number of BH 
detections has increased to the degree that the field has 
shifted from confirming the existence of BHs to compar- 
ing their properties with those of the host galaxies. One 
of the most discussed relationships found between BHs 
and host galaxies is the correlation between BH mass and 
the velocity dispersion of stars in the bulge (ie, the Mbh- 
(7* relation; Ferrarese & Merritt (2000); Gebhardt et al. 
(2000); Giitekin et al. (2009), and references therein). 
It is often suggested that this relation is a manifestation 
of a causal connection between the formation and evo- 
lut ion of the black hole a nd the bulge. As pointed out 
bv iGebhardt et ail (|2000h . it is natural to assume that 
bulges, black holes, and quasars formed, grew, or "turned 
on" as parts of the same process, in part because the col- 
lapse or merger of bulges might provide a rich fuel supply 
to a centrally located black hole. However, the nature of 
this connection is still not well understood. 

Clearly, accurate BH mass measurements in mega- 
maser galaxies can help to improve our understanding of 
the Mbh-o* correlation, especially at the low-mass end. 
We present high-resolution images and rotation curves 
of megamaser disks in six new galaxies and their BH 



masses. In addition, we r e-analyze the BH m ass in UGC 
3789 using the data from lReid et al.l (|2009aD and present 
the new measurement here. In section 2, we present our 
sample of galaxies, VLBI observations and data reduc- 
tion. Section 3 shows the VLBI images and rotation 
curves of the megamaser disks, followed by our analysis 
of BH masses. In section 4, we rule out compact clus- 
ters of stars or stellar remnants as the central dominant 
masses for the majority of the megamaser disks. In sec- 
tion 5, we compare results from other mass measuring 
techniques to our maser masses. The discussion of our 
main results is presented in section 6, and we summarize 
the results in section 7. A more in-depth discussion of the 
Mbh-c* relation, including the ne w maser BH masses, i s 
presented in a companion paper bv I Greene "eFaD (l2010h . 

2. THE SAMPLE, OBSERVATION, AND DATA REDUCTION 

2.1. The Megamaser Disk Sample 

Water masers have so far been detected in ~ 1 36 galax - 
ies Q Most of the megamasers originate in AGN (|Loll2005l) 
and about 20% of them show spectra suggestive of emis- 
sion from sub-parsec scale, edge-on, circumnuclear disks. 
In the MCP, we are currently conducting VLBI obser- 
vations of six megamaser disks to determine the angular 
diameter distances to their host galaxies. While more 
data are needed to measure their distances accurately, 
the quality of the VLBI imaging is already good enough 
to measure precise central BH masses. In addition to 
the six galaxies for distance measurements, we also have 
VLBI data on the megamaser disk in NGC 4388 to mea- 
sure the central BH mass, even though the data are not 
suitable for a distance determination. Table 1 lists coor- 
dinates, recession velocities, spectral types, and morpho- 
logical types for these seven galaxies. 

2.2. Observations 

The megamaser galaxies in our sample were observed 
between 2005 and 2009 with the Very Long Baseline 
Array (VLBA)0, augmented by the 100-m Green Bank 
Telescope (GBT) and in most cases the Effelsberg 100-m 
telescope^ Table 2 shows the basic observing informa- 
tion including experiment code, date observed, antennas 
used, and sensitivity. 

We observed the megamasers either in a phase- 
referencing or self-calibration mode. With phase- 
referencing we perform rapid switching of the telescope 
pointing between the target source and a nearby (< 1°) 
phase calibrator (every ~ 50 seconds) to correct phase 
variations caused by the atmosphere. In a self-calibration 
observation, we use the brightest maser line(s) to cali- 
brate the atmospheric phase. In both types of observa- 
tions, we placed "geodetic" blocks at the beginning and 
end of the observations to solve f or atmosphere and clock 
delay residuals for each antenna (jReid et al.ll2009bh . For 
NGC 6323 and NGC 6264, we also placed a geodetic 
block in the middle of the observations to avoid the 

1 https: / /safe. nrao.edu/wiki/bin/view/Main/MegamaserCosmologyProjcct 

2 The VLBA is a facility of the National Radio Astronomy Ob- 
servatory, which is operated by the Associated Universities, Inc. 
under a cooperative agreement with the National Science Founda- 
tion (NSF). 

3 The Effelsberg 100-m telescope is a facility of the Max-Planck- 
Institut fur Radioastronomie 
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zenith transit problem at the GBT and to obtain bet- 
ter calibration. In each geodetic block, we observed 12 
to 15 compact radio quasars that cover a wide range of 
zenith angles, and we measured the antenna zenith de- 
lay residuals to ~ 1 cm accuracy. These geodetic data 
were taken in left circular polarization with eight 16-MHz 
bands that spanned ~ 370 - 490 MHz bandwidth centered 
at a frequency around 22 GHz; the bands were spaced 
in a "minimum redundancy" manner to sample, as uni- 
formly as possible, all frequency differences between IF 
bands in order to minimize ambigui ty in the delay so- 
lution (jMioduszewski fc Kogari l2009\ We also observed 
strong compact radio quasars every ~ 20 minutes to 2 
hours in order to monitor the single-band delays and elec- 
tronic phase differences among and across the IF bands. 
The errors of the single-band delays are < 1 nanosecond. 

2.3. Data Reduction 

We calibrated all the data using the NRAO As- 
tronomical Image Processing System (AIPS). Since 
the geodetic data and the maser data used different 
frequency settings, we reduced them separately. For 
the geodetic data, we first calibrated the ionospheric 
delays using total electr on content measurements 
(|Walker fc Chatterieell2000f) and the Earth Orientation 
Parameters (EOPs) in the VLBA correlators with 
the EOP estimates from the US Naval Observatory 
(|http:/ /gemini.gsfc.nasa.gov/solve-Save/usnoJinals.erp[ ). 
We performed "fringe fitting" to determine the phases, 
single-band delays, and fringe rates of IF bands of each 
antenna for every geodetic source, and the multi-band 
delay of each antenna was determined from these 
solutions. Finally, we determined the residual tropo- 
spheric delay and clock errors for all antennas using the 
multi-band delays. We applied these corrections to the 
maser data as described in the next paragraph. In five 
datasets for NGC 6323, we also made antenna position 
corrections before the fringe-fitting process when the 
positions used in the observation deviated from the 
latest USNO (http://rorf.usno.navy.mil/solutions/) 
solutions by more than ~ 2 cm. 

For the maser data, after the initial editing of bad data, 
we corrected for ionospheric delay and the EOPs in the 
same way as for the geodetic dataset. We then corrected 
the sampler bias in the 2-bit correlator. The amplitude 
calibration was done with the information in the gain ta- 
ble and the system temperature table. We corrected the 
interferometer delays and phases caused by the effects of 
diurnal feed rotation (parallactic angle), and applied the 
tropospheric delay and clock corrections obtained from 
the geodetic data afterwards. The next step was to per- 
form fringe fitting on one or two scans of the delay cal- 
ibrators to calibrate the electronic phase offsets among 
and across IF bands. The frequency axes of the maser in- 
terferometer spectra were then shifted to compensate for 
the changes in source Doppler shifts over the observing 
tracks. 

The final step in calibration was to solve for the atmo- 
spheric phase variation by using either phase-referencing 
or self-calibration. In the phase-referencing mode, we 
ran the AIPS task CALIB on the phase calibrator to de- 
termine the phase correction as a function of time for 
each individual IF band. In the self-calibration mode, 
we either selected a single strong maser line or averaged 



multiple maser lines in narrow ranges of both velocity 
and space to perform phase calibration. The typical so- 
lution interval was 100 seconds. After the above calibra- 
tions, we discarded the phase solutions and the maser 
data in the time intervals within which the solutions ap- 
peared to be randomly scattered in time or had adja- 
cent reference phases exceeding ~50°. The phase solu- 
tions were then interpolated and applied to all the maser 
data. In all phase-referencing observations for NGC 
6323, we performed an extra step of self-calibration af- 
ter the initial phase-referencing calibration, because the 
phase-referencing calibration alone did not give adequate 
phase calibration, and significant sidelobes were still vis- 
ible after CLEAN deconvolution. 

The calibrated data from multiple tracks of a single 
source were combined for imaging. We Fourier trans- 
formed the gridded (u, v) data to make images of the 
masers in all spectral channels of the IF bands that 
showed maser lines, and we deconvolved the images using 
CLEAN with a weighting scheme that optimized the po- 
sition accuracy (e.g. ROBUST = in the task IMAGR). 
We fitted the detected maser spots with elliptical Gaus- 
sians to obtain the positions and flux densities of individ- 
ual maser components. In Table 3 we show a represen- 
tative dataset that includes the velocities, positions, and 
peak intensities of the maser spots in NGC 6264. The 
data for all galaxies are available in the electronic version 
of this paper. Note that all galaxies except NGC 4388 
and NGC 2273 have at least 2 tracks of data. There- 
fore, the peak flux in Table 3 is the flux for the averaged 
dataset. The actual fluxes in an individual epoch can be 
higher or lower because of variability of the masers. 

2.4. Relativistic Velocity Assignment 

The sub-parsec megamaser disks presented in this pa- 
per are in a deep gravitational potential and the majority 
have recession velocities over 1% of the speed of light c. 
For this reason, we made both special and general rel- 
ativistic corrections to the maser velocities in the data 
before we used the data to analyze the BH masses. We 
made the relativistic corrections in the following way: 

The observed Local Standard of Rest velocities V op 
listed in Table 3 are based on the "optical" velocity con- 
vention 

where v is the observed frequency and vq = 22.23508 
GHz, the rest frequency of the H2O 616 — 523 transition. 

Because of gravitational time dilation, the emitting fre- 
quency vq of a maser at distance r from a compact object 
of mass M and the actual observed frequency z/qo for an 
observer at r = 00 differ by a factor 
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rc z I 



(2) 



For a maser in a circular orbit moving at speed /3 m c, 
balancing the gravitational and centripetal accelerations 
gives GM/r 2 — (/3 m c) 2 /r, so 



— = 1 + P 2 

is, + Pn 



(3) 



We multiplied the observed frequency v of each maser 
line by (l+/3^) to correct for gravitational time dilation. 
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The megamaser lines at the systemic velocities were 
also corrected for a transverse Doppler shift 



X/2 



(4) 



The relativistically correct Doppler equation for a 
source moving radially away with velocity v — /3c is 



1-/3 
1 + /3 



1/2 



After the observed frequency had been corrected for 
the gravitational time dilation and transverse Doppler 
shift, we used Equation [5] to convert the corrected fre- 
quency to its relativistically correct radial velocity v = 
f3c. For the megamaser disks, typically (3 m < 0.003. So, 
for galaxies in the Hubble flow (z > 0.01), the general 
relativistic corrections are smaller than the special rel- 
ativistic corrections. For example, for masers in UGC 
3789, which has an optical-LSR recession velocity of 3262 
km s -1 , the special relativistic corrections range from 10 
to 26 km s _1 whereas the general relativistic corrections 
range from to ~ 2 km s _1 . 

Finally, when fitting the rotation curves, we used the 
relativistic formula for the addition of velocities to de- 
compose the observed j3 value of each maser spot into 
a common /3 g associated with the radial velocity of the 
center of each megamaser disk and an individual j3 m as- 
sociated with orbital motion of a specific maser spot. 
The redshifted maser spots have (3 m > and blueshifted 
maser spots have f3 m < 0. 



= 



1 + A^n 



3. RESULTS 



3.1. VLBI Images, Rotation Curves, and BH Masses 

Figure 1 shows the GBT single-dish spectra for all 
megamaser galaxies prese nted here ex c ept for UGC 3789, 
which can be found in IReid et al.l (|2009aj ). Figures 
2, 3, and 4 show the VLBI maps and the position- 
velocity (P — V) diagrams along with the fitted rota- 
tion curves of the maser spectral components (spots) 
in UGC 3789, NGC 1194, NGC 2273, NGC 2960 
(Mrk 1419 ) NGC 4388, NGC 6264, and NGC 6323. 
IReid et al.l (|2009al) published the UGC 3789 VLBI map 
and the rotation curve, and we performed a new analysis 
of the BH mass for this galaxy based on those data. We 
show the VLBI map and rotation curve for this galaxy 
again for direct comparison with the other six megamaser 
disks. The data points in the VLBI maps and rotation 
curves are color-coded to indicate redshifted, blueshifted, 
and systemic masers, where the "systemic" masers refer 
to the maser spectral components having velocities close 
to the systemic velocity of the galaxy. Except for NGC 
4388, the maser spot distributions are plotted relative to 
the average position of the systemic masers. Systemic 
masers are not detected in NGC 4388, so we plotted its 
maser distribution relative to the dynamical center de- 
termined from fitting the data in the P—V diagram with 
a Keplerian rotation curve. 

To estimate the inclination and dynamical center (i.e. 
the position of the BH) of each disk, we rotated the coor- 
dinate system to make the disk horizontal and used the 



fitted horizontal line that passes through the high veloc- 
ity masers as the zero point of the y-coordinate of the 
dynamical center (see Figures 3 & 4). The zero point of 
the x-coordinate is defined to be the unweighted average 
8 X of the systemic masers. 

We assumed that the systemic masers have about the 
same orbital radii as the high-velocity masers to esti- 
mate the maser disk inclination cos 

~\{4" V ' ) )/ 6 t), where 

(5) (@y yS ^) is the average 9 y position of the systemic masers 
and 8 r is the orbital radius of the systemic masers. 
In principle, one can determine # r exactly only when 
good rotation curves for both systemic and high- velocity 
masers can be obtained, and precise acceleration mea- 
surements for the systemic masers are available. Among 
our data we only have such info rmation for UGC 3789 
at this point (jBraatz et al.ll2010h . and so we use (9 X ) of 
high velocity masers as an estimate of r . We note that 
all our megamaser disks with systemic masers detected 
have inclinations larger than 80°, so assuming the disk is 
exactly edge-on will only cause errors less than 1% in the 
derived BH masses. In NGC 4388 we could not measure 
a disk inclination, but even if the disk were 20° from 
edge-on, the contribution to the error in the BH mass 
would be only 12%, comparable to the 11% error caused 
by the distance uncertainty. 

We determined the rotation curve for each megamaser 
disk as a function of the "impact parameter" defined as 
the projected radial offset 9 = (8 2 + O 2 ) 1 / 2 of the maser 
spots so that we can account for the warped structures 
in some megamaser disks. We then performed a nonlin- 
ear least-squares fit of a Keplerian rotation curve to the 
position-velocity diagram with the assumption that the 
high-velocity masers lie exactly on the mid-line of the 
disk. In addition, the systemic velocity of each galaxy 
was fitted as a free parameter, and we report the best 
fits of the systemic velocities of our megamaser galaxies 
in Table 1. 

The fitted Keplerian rotation curves can be expressed 
with the following form: 



(6) 



\VK \ = Vi 



1 mas 



-1/2 



(?) 



where |uk| is the orbital velocity (after relativistic cor- 
rections) of the high velocity masers and v\ is the orbital 
velocity at a radius 1 mas from the dynamical center. 
The BH mass is 



M. 



G 



D A 



6.48 x 10 8 G 



(8) 



where Da is the angular diameter distance to the galaxy. 
We show all the measured BH masses in Table 4. 

3.2. The Error Budget For the BH Mass 

There are three major sources of error for our BH 
mass calculation, and the largest comes from the dis- 
tance uncertainty. Except for NGC 4388, we used the 
Hubble distances relative to the CMB for a cosmological 
model with H = 73 km s" 1 Mpc -1 , f2 ma ttcr = 0.27 and 
17a = 0.73. Since Hp is uncertain by 8 km s _1 Mpc -1 
(jFreedman et al.l 120011) . the galaxy distances will have 
an rms uncertainty of 11%. NGC 4388 is in the Virgo 
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cluster and has a large peculiar velocity, so w e used its 
Tully-Fisher distance of 1 9 Mpc ([ Russell 2002). Since the 
distance modulus used bv lRusselll ([20021 ) is calibrated rel- 
ative to the Cepheid distances prov ided by the Hubble 
Key Project (jFreedman et al.ll2001| ) and the systematic 
error of the Cepheid distances is the primary source of 
the Hubble constant uncertainty, we also adopted 11% 
as the rms distance error for NGC 4388. 

The second source of BH mass uncertainty comes from 
fitting the rotation curves of the masers. In our Keple- 
rian fitting, the rms o~g of the observed position offsets 
from Keplerian rotation is usually larger than the rms 
position uncertainties of the data by a factor of 1.5 to 4. 
This extra scatter can lead to a systematically different 
position of the dynamical center, and hence a slightly 
different BH mass if we allow dynamical center position 
to be a free parameter in the fit. This excess may in- 
dicate that we underestimated our observational errors 
or it may indicate genuine deviations from our simple 
model, possibly caused by peculiar motions of the masers 
(e.g. perturbations from spiral density waves) or high- 
veloc ity masers not lying prec isely on the mid-line of the 
disk ([Humphreys et al.l 120081 ). If the latter is the mam 
cause of deviation, we estimate that the majority of the 
maser spots fall within ss 7° — 10° from the midline of 
the disk, depending on the galaxy. 

In addition to the deviations of the masers from the 
midline of the maser disk, there is also an error in the 
BH mass from the uncertainty of the position of the BH 
in the fitting. To estimate this error, we relax our as- 
sumption on the position of the BH, and allow it to be 
a free parameter in the fitting. However, since the reces- 
sion velocity of the galaxy (v ) and the position of the 
BH are correlated in this case, we impose a constraint 
on the possible positions of the BH such that the fitted 
vo does not exceed vo determined from other meothods 
(e.g. HI measurements) beyond their error bars. Includ- 
ing all possible errors mentioned above, we estimate the 
total fitting error in BH mass to be 1%— 5% depending 
on the galaxy. 

The third source of error is from the absolute posi- 
tion errors of our megamasers. Source position errors 
introduce an extra phase difference between maser spots 
having different frequencies (velocities) . The relative po- 
sition errors sc ale with velocity offset from the referenc e 
maser feature ([Argon et alJl2007t iGreenhill et al. 1993), 
which may differ from source to source. Since the BH 
mass is proportional to the size of the maser orbits, rel- 
ative position shifts among masers introduce errors in 
the BH mass measurements. Among our megamasers, 
UGC 3789, NGC 2960, and NGC 2273 have the largest 
absolute position errors (10 mas; see Note (h) in Table 
1). The resultant BH mass errors are 3 — 5%, whereas 
this error is smaller than 0.3% in the other megamaser 
galaxies. Please note that the error for the BH mass dis- 
cussed in the following subsections does not include the 
distance uncertainty. Only source position errors and er- 
rors from Keplerian fitting are evaluated, which is the 
best approach for comparing our results with BH masses 
derived from other techniques as the same distance can 
be used for comparison. As the distance uncertainties 
shrink in the future, so will the uncertainties in our BH 
masses. 



3.3. Notes On Individual Galaxies 

3.3.1. NGC 1194 

NGC 1194 hosts a Seyfert 1.9 nucleus and has a dis- 
tance of « 52 Mpc. The position angle of the maser disk 
(Figure 2) is 157° and the inclination is sa 85°. NGC 
1194 has the largest maser disk among the megamasers 
presented here, with an inner and outer radius of 0.53 
and 1.3 pc, respectively. The blueshifted and redshifted 
masers do not appear to fall exactly on a straight line 
on the sky, but they are consistent with a slightly bent 
thin disk. The rotation curve (Figure 3) is consistent 
with Keplerian rotation, and the measured BH mass is 
(6.5 ± 0.3) x 10 7 M Q , the largest among the BH masses 
studied here. 

3.3.2. NGC 2273 

NGC 2273 is a Seyfert 2 galaxy at a distance of « 26 
Mpc. The position angle of the maser disk (Figure 2) is 
153° with an inclination of « 84°. The inner and outer 
radii of the disk are 0.09 and 0.26 pc, respectively. The 
maser distribution shows a hint of a warped disk. The ro- 
tation curve (Figure 3) is also consistent with Keplerian 
rotation to within the errors. The BH mass we obtain is 
(7.6 ±0.4) x 10 6 M Q . 

3.3.3. UGC 3789 

UG C 3789 is a Seyfer t 2 galaxy at a dista nce of « 50 
Mpc dReid et al.l l2009at iBraatz et al.ll20loh . Its edge- 
on maser disk lies at position angle 41° East of North 
(Fig. 2). Although the high- velocity masers appear to 
trace a flat disk, the misalignment of the systemic masers 
with respect to the disk and the fact that the systemic 
maser s have at least tw o distinct centripetal accelera- 
tions ([Braatz e t al. 2010) suggest that the disk is proba- 
bly warped along the line of sight. The inner and outer 
radii of the disk are 0.085 and 0.32 pc, respectively. Fit- 
ting the P—V diagram (Fig. 3) with a Keplerian rotation 
curve gives a BH mass of (1.12 ± 0.05) x 10 7 A/ Q . Note 
that this ma s s is n early the same as the mass given in 
IBraatz et all (pOlOl ). (1.09 ± 0.05) x 10 7 M Q . The very 
small difference comes from choosing different positions 
as the dynamical center, the uncertainty of which has 
been taken into account in our estimate of the BH mass 
error. 

3.3.4. NGC 2960 (Mrk 1419) 

NGC 2960 hosts a LINER nucleus at a distance of 71 
Mpc. The position angle of the maser disk (Fig. 2) is 
— 131° and the inclination is ~ 89°. The outer parts 
of the disk show some warping. The inner and outer 
radii of the disk are 0.13 and 0.55 pc, respectively. The 
rms scatter of the high velocity masers normal to the 
disk is o~± = 71/iarcsec, which is 1.6 times larger than 
the rms uncertainty of the data, so it is likely that ei- 
ther we underestimated the observational uncertainties 
because of larger tropospheric delay errors for low decli- 
nation sources, or the thickness of the disk may not be 
negligible in this megamaser. In this galaxy, we measure 
a BH mass of (1.14 ± 0.05) x 10 7 M Q . 

3.3.5. NGC 4388 

NGC 4388 is in the Virgo cluster and we adopted a 
Tully-Fisher distance of 19 Mpc; it is the nearest of the 
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seven galaxies presented here. It hosts a Seyfert 2 nu- 
cleus and the maser disk (Figure 2) has a position angle 
of 107°. We cannot measure the inclination directly be- 
cause no systemic masers were detected. As with the 
other megamasers, we assumed a flat, edge-on disk. 

Since no systemic masers were detected, we determined 
the position of the dynamical center by allowing it to be 
a free parameter when fitting the high-velocity masers 
with a Keplerian rotation curve (Fig. 4). In addition, 
we fixed the s ystemic velocity of the galaxy using an HI 
measurement (jLu et al.ll2003t ) . The resulting BH mass is 
(8.5±0.2) x 10 6 M Q . With only five maser spots mapped, 
there is not sufficient data to demonstrate Keplerian ro- 
tation or even to show that the masers actually lie on a 
disk. We argue that it is l ikely they do because the radio 
continuum jet shown by Kuku laet al.l (|l995) is nearly 
perpendicular to the line joining the blueshifted and red- 
shifted masers. However, because of these uncertainties, 
the BH mass for NGC 4388 should be used with some 
caution until better data are obtained. 

3.3.6. NGC 6264 

At w 136 Mpc, NGC 6264 (Figure 2) is the most dis- 
tant object in our current sample. The disk is slightly 
warped and appears to have some thickness, but given 
the uncertainties in the position measurement, our ob- 
servation is consistent with a thin disk. The disk has a 
position angle of —85° and an inclination of ss 90°. The 
inner and outer radii of the disk are 0.18 and 0.77 pc, 
respectively. In the P — V diagram (Fig. 4), the high 
velocity masers beautifully trace the Keplerian rotation 
curve, and we obtain a BH mass of (2.84±0.04) x 10 7 M Q . 

3.3.7. NGC 6323 

NGC 6323 is a Seyfert 2 galaxy at a distance of « 105 
Mpc. The VLBI image (Fig. 2) shows a remarkably thin 
disk at a position angle of 10° and inclination of s» 89°. 
The disk is apparently warped, and the inner and outer 
radii are 0.13 and 0.3 pc, respectively. The rotation curve 
of the high- velocity masers (Figure 4) is Keplerian, and 
we obtain a BH mass of (9.3 ± 0.1) x 1O 6 M . 

3.4. Search For Continuum Emission 

We searched for continuum emission from the vicinity 
of the supermassive BH (i.e., near the systemic masers) 
in our megamaser galaxies by averaging the line-free 
spectral channels in our data and imaging with natu- 
ral weighting to maximize the detection sensitivity. We 
detected no continuum emission in all megamaser galax- 
ies presented in this paper. The channels averaged, the 
center velocities of the bands used for averaging, and the 
continuum upper limits are listed in Table 5. 

4. A SUPERMASSIVE BLACK HOLE OR A CENTRAL 
CLUSTER OF STARS OR STELLAR REMNANTS ? 

The near-Keplerian shapes of the rotation curves of 
many circumnuclear megamaser disks allows one to de- 
termine the enclosed (spherical) mass within the inner- 
most radius of the disk. While the enclosed masses of 
the megamaser disks are very likely dominated by super- 
massive BHs, as we assumed in the previous sections, it 
is still important to find ways to justify this assumption. 
The question of whether the enclosed mass can be at- 
tributed to a point mass or a compact cluster of stars or 



stellar remn ants has been addressed by iMaozl (|1995| ) & 
iMaozl (|1998l ). The main argument is that if the lifetime 
of a central cluster, limited by evaporation or collision 
timescales, is significantly shorter than the age of its host 
galaxy, then it is unlikely to persist, and a central super- 
massive BH would be required to account for the enclosed 
mass. Here, the evaporation timescale t. 



cvap 



136* 



relax 



(^roiax is the half-mass relaxation timescale; Binney & 
Tremaine 2008) is considered to be the upper limit of 
the lifetime of any bound stellar system whereas the col- 
lision timescale is the characteristic timescale that a star 
suffers a physical collision (i.e. an inelastic encounter; 
Binney & Tremaine 2008). 

To estima te the lifetime of the central cluster, we follow 
iMaol ((1991 and assume that the central cluster has the 
Plummer density distribution (|Plummerlll915D : 



p(r) = Po 1 + 



-5/2 



(9) 



where po is the central density and r c is the core radius. 
(The reason for choosin g the P l umme r distribution is de- 
scribed in section 2.1 of lMaozl (|1998|) ) 

We constrain po and r c by fitting the position- velocity 
diagram of the megamaser disk with the rotation curve 
of a Plummer cluster: 



Vp 



GM a 



_|_ 
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(10) 



where — Anporl/S is the total mass of the cluster. 
Here, M x and r c are fitted as free parameters. From 
these two parameters, we calculated po = 3M 00 /47rr^. 
In all cases is very close to the "enclosed" mass 
measured from the Keplerian rotation curve fit in Section 
3, and the differences are less than 4 — 18%. In Figures 3 
& 4, we show the fitted Plummer rotation curves along 
with the fitted Keplerian rotation curves. 

Rather than using the method described above to con- 
strain the core radius r c , in some cases we could apply 
another approach that places even tighter constraints. 
We note that the Plummer rotation curve does not de- 
crease monotonically with radius; instead the rotation 
curve turns over at a maximum rotation speed v max : 



( 2GM 



V 3 3 / 2 



1/2 



(11) 



Having a maximum rotation velocity is not unique to 
a Plummer cluster. It is a general feature for clusters 
having the same form of density profile with the exponent 
smaller than —3/2. 

The core radius r c of a Plummer sphere having maxi- 
mum rotation speed v max 1S 



.(max) 



2 GM 00 

2 2 
max 



(12) 



We used Equation 12 to estimate r c for the megamaser 
disks in NGC 1194, NGC 2273, NGC 2960 and NGC 
4388. In these cases, we do not have very well-sampled 
or high quality rotation curves (see Figs. 3 & 4), so using 
r (max) £ r0ln Eq ua tion 12 actually sets a tighter constraint 
on the core radius of the cluster than using rotation-curve 
fitting. For these four cases, we use the highest observed 



7 



velocity in the PV diagram as an estimate (lower limit) 

of v max and use it to calculate rc max ^ with from the 
Plummer rotation curve fitting from Equation 10. For 
UGC 3789, NGC 6323, and NGC 6264, we used r c from 
rotation curve fitting because the qualities of the rotation 
curves are good, and they give tighter constraints on r c . 
In Table 4, we give the Plummer model parameters for 
all of our megamaser disks. 

We constrained the lifetime of the Plummer cluster 
T age by first requiring the cluster not evaporate in a 
timescale less than the age of its host galaxy (> 10 Gyr 
if the galaxy has formed before z = 2). This requirement 
sets an upper limit to the mass of the constituent stars of 
the cluster because evaporation is unimportant so long 
as the mass of its stars satisfies the following equation: 



TO* 



< 



KAMoo/m*)]- 



3/2 



( M 



1/2 



m q j ~ / —/j yo.oipcy V 108M © 

'(13) 

where to* is the mass of each star, A pa 0.1 
(|Binnev fc Trem ainc 20QH), and is the radius of half 
total mass (r^ pa 1.305r c ). We call the maximum m* that 
satisfies the above equation m max (Table 4) and used it to 
calculate the collision timescale of the Plummer cluster 
with the Plummer model parameters: 



icon 
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where n* = po/m max is the number density of stars, c* 
is the rms vel ocity dispersion of the stars, and r* is the 
stellar radius (jMaozlll998l; iBinnev fc Tremaind 120081 ) . If 
tcoll < 10 Gyr, then T age is constrained by t co u and we 
can rule out the Plummer cluster as an alternative to 
the BH. If t co u > 10 Gyr, we cannot rule out a cluster 
whether r age is dominated by evaporation or collision. 

For UGC 3789 we obtained an upper limit m max ~ 
0.12M© to the mass of individual stars and a lower limit 
N > 1.0 x 10 8 to the number of stars in the cluster. 
The mass limit directly rules out neutron stars as the 
constituents of the cluster, and only brown dwarfs, very- 
low mass stars, or white dwarfs are possible. If the 
constituents of the cluster are brown dwarfs, the colli- 
sion timescale for our Plummer model in UGC 3789 is 
t C o\\ < 5.0 x 10 6 years, much less than the age of a galaxy. 
The timescale is even shorter if the constituent stars are 
main-sequence stars. However, if the cluster is composed 
of white dwarfs, the collision timescale can be as long as 
2.7 x 10 7 years, but this is still much shorter than the life- 
time of a galaxy. Therefore, we conclude that a compact 
cluster is not likely to survive long, and the dominant 
mass at the center of UGC 3789 is a supermassive BH. 
For all other megamaser galaxies except NGC 1194, the 
lifetimes of the central clusters are also shsorter than 
the age of a galaxy (Table 4). The constraints on life- 
time weakly rule out massive clusters in NGC 4388 and 
NGC 6264 ( T agc w 7 x 10 9 yr) but strongly rule out 
clusters in the others (r age < 10 9 yr). 

In summary, by setting tight constraints on the sizes 
and central mass densities of possible Plummer clus- 
ters, we have been able to strongly rule out clusters as 
the dominant central masses in UGC 3079, NGC 2273, 
NGC 6323, and NGC 2960 and weakly rule out the clus- 



ters in NGC 4388 and NGC 6264. We argue that 
supermassive black holes are the dominant masses in 
these megamasers. To gether with the M ilky Way Galaxy, 
NGC 4258, and M3 1 |Korme"ndvl [200l . the number of 
galaxies with strong evidence to rule out a massive star 
cluster as the dominant central mass increases from three 
to seven. 

5. OTHER BH MASS MEASURING TECHNIQUES 

It is important to compare results from different meth- 
ods for measuring BH masses, because such compar- 
isons can provide in sight into potential systematic er- 
rors for each method (ISiopis et al.ll2009t IKormendvll2004l ; 
Humphrey etail 120091: IGreene et al.l |2010TT Compar- 
ing optically measured dynamical BH masses with those 
from the H2O megamaser method is especially valuable, 
since the megamaser galaxies with Keplerian rotation 
curves provide the most direct and accurate BH mass 
measurements for external galaxies, These maser BH 
masses can be used to test the more commonly used BH 
mass measuring techniques in the optical, such as the 
stellar or gas dynamical methods. 

One cannot meaningfully compare the stellar or gas dy- 
namical BH masses with our maser BH masses unless the 
gravitational spheres of influence can be resolved. Other- 
wise, even if the optically determined dynamical masses 
agree with the maser BH masses within the errors, the 
uncertainty can be too large to tell the real accuracy of 
the stellar- or gas-dynamical method. The angular diam- 
eters of the spheres of influence of our maser BH masses, 
2.Ri nf (Table 4), range from 0'.'006 to 0'.'06. Among these 
megamasers, only the spheres of influence in NGC 4388 
and NGC 1194 can be barely resolved by the HST (res- 
olution « 0'.'07 at A sa 65001) or the VLT assisted 
with adaptive optics (resolution pa O'.'l at A = 2/im). 
Therefore, of the galaxies here, stellar or gas dynamical 
measurements are feasible only for these two galaxies. 
We have obtained VLT time to measure the BH mass 
in NGC 4388, and we will apply the stellar dynamical 
method to this galaxy and compare the BH mass to the 
maser BH mass in the future. 

Another commonly used BH ma ss measuring tech- 
nique is the viria l estimation method (IGreene 
iKim et al.l 120081: iVestergaard fc Osmerl 120091) . We are 
able to compare this technique with the megamaser 
disk method in four galaxies, and this is the first time 
that the virial estimation method is directly tested. In 
this method, one estimates the BH mass as Mbh = 
/ RBLRV? ine /G, where / is an unknown factor that de- 
pends on the structure, kinematics, and orientation of 
the broad-line region (BLR), Rblr is the radius of the 
BLR, (Jii ne is the gas velocity dispersion observed in the 
BLR, and G is the gravitational constant. Here, we adopt 
the latest e mpiri c ally determin e d (/) = 5.2 ± 1.3 from 
IWoo et all 1)20101) . IVestergaard] (|2009f) suggests that the 
virial method is accurate to a factor of pa 4. Since we 
cannot directly detect the BLRs in megamaser galaxies 
that have Seyfert 2 nuclei, we estimated au ne using the 
scattered "polarized broad lines" (PBL) from the hidden 
BLRs in those four megamaser galaxies with detected 
PBLs: NGC 1068, NGC 4388, NGC 2273, and Circinus. 
We estimated -Rblr via the L( 2 -iokev) —Rblr correlation 

4 The Very Large Telescope 
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(|Kaspi et a l. 2005). Table 6 shows the resulting virial BH 
masses for these four galaxies. 

Under the assumption that the observed linewidth of 
broad H a or Ha emission approximates the "intrinsic 
linewidth'fl, we find that for NGC 1068, NGC 2273, 
NGC 4388 and Circinus, the BH masses at the 1 a level 
measured by the virial estimation method agree within 
a factor of 5 with the megamascr BH masses. So, the 
scatter of the virial BH masses relative to the maser BH 
masses is comparable to the factor of 4 accuracy expected 
for the virial method. It is possible that the slightly 
larger difference between the maser and viral BH masses 
as compared with the expected accuracy for the virial 
method is partly caused by the larger intrinsic scatter 
in the £(2-iofceu) — ^?blr correlation which we used here 
than that in the more generally used \L\ (5100A) — Rblr 
correlation (e.g. Kaspi et al. 2005) in the optical. 

The application of the virial method to megamaser 
galaxies is described in detail in the Appendix. 

6. MASER BH MASSES AND THE M - a* RELATION 

Table 4 shows the BH masses for all seven megamaser 
galaxies. Along with the previously published BH masses 
from megamaser observations, it is interesting to see that 
the maser BH masses, except in Circinus, are within a 
factor of 3 of 2.2 x 10 7 M©, with the uncertainty of each 
BH mass < 12%. This small range of masses could be 
a selection bias. Disk megamasers are preferentially de- 
tected in Seyfert 2 spiral galaxies, and since the local 
active BH mass funct ion for Seyfert 2 gala xies peaks at 
M BH ~ 3 x 10 7 M Q (|Heckman et al.ll200l . we may be 
limited to BHs in this range of masses by analyzing disk 
megamasers. 

Our new maser BH masses more than double the num- 
ber of galaxies having dynamical BH masses A^bh ~ 
10 7 Mq. These measurements play a particularly impor- 
tant role in constraining the M — er* relation at the low- 
mass e nd of known nuclear BH masses. In a companion 
paper, iGreene et al.l (|2010l ) find that the maser galaxies 
as a group fall significantly below the M — a* relation de- 
fined by more massive elliptical galaxies. As a result, the 
Al — (j* relation that fits later- type and lower-mass galax- 
ies has both a larger scatter and lower zeropoint than the 
relation for elliptical galaxies alone. However, there is a 
potential caveat that measuring rob ust (well-defined) a y 
in late-type galaxies is challenging ([Greene eTaD[2Qi3), 
and the contribution of systematic biases in er* to the de- 
viation from the M—a± relation still needs to be explored 
in the future. With this caveat in mind, the observed 
deviations from the M — a* relation at low mass imply 
that the relation may not be a single, low-scatter power 
law as originally proposed, and our BH mass distribution 
is consistent with the best simulations of galaxy growth 
regulated by "radio" feedback; see Figure 4 of Croton et 
al. (2006). 

In addition to our seven galaxies with maser BH 
masses, we currently have VLBI datasets for another four 
disk megamaser candidates, and we are monitoring the 
spectra of more than six disk megamaser candidates that 
are currently too faint to be observed with VLBI, but 
may flare up in the future. Along with these galaxies and 

5 By "intrinsic linewidth", we mean the linewidth one would 
measure if the BLRs were observed directly as in Type 1 AGNs 



more megamasers we may discover in the future with the 
GBT, we expect to obtain a larger sample of maser BH 
masses which could help to clarify the A/bh-c* relation. 

7. SUMMARY 
Our main conclusions are the following: 

1. The maser distributions in all seven megamaser 
galaxies are consistent with edge-on circumnuclear 
disks surrounding central massive objects in the ac- 
tive galactic nuclei. The inner radii of the disks are 
between 0.09 and 0.5 pc, similar to all previously 
published megamaser disks. The rotation curves of 
all seven megamaser disks are consistent with Kep- 
lerian rotation. Four of the megamaser disks reveal 
evidence for warps. 

2. VLBI observations of circumnuclear megamaser 
disks are the only means to measure directly the 
enclosed mass and the mass density well within the 
radius of the gravitational sphere of influence of the 
central mass in these galaxies. The high central 
mass densities (0.12 to 60 xlO lo M pc~^) (within 
the central 0.3 pc) of the seven megamaser disks 
indicate that in all, except one, disks, the central 
mass is dominated by a supermassive BH rather 
than an extremely dense cluster of stars or stellar 
remnants. 

3. The BH masses measured are all within a factor of 
3 of 2.2 x 10 7 A/© and the accuracy of each BH mass 
is primarily limited by the accuracy of the Hubble 
constant, currently ±11%. The narrow range of 
BH mass distribution may reflect selection from the 
local active-galaxy mass function. 

4. Stellar dynamics cannot be applied to the major- 
ity of the megamasers presented here to measure 
BH masses with high precision. The gravitational 
spheres of influence in all cases except NGC 1194 
and NGC 4388 are too small to be resolved by 
current optical telescopes. Observations with the 
VLT for measuring the BH mass in NGC 4388 with 
stellar-dynamical modeling are in progress. 

5. Under the assumption that the broad emission 
linewidths can be estimated from polarized scat- 
tered light, we have calibrated for the first time 
the BH mass determination by the virial estima- 
tion method based on optical observations. With 
the latest empirically determined (/) = 5.2 ± 1.3, 
the virial estimated BH mass is within a factor of 5 
at the 1 a level of the accurate BH mass based on 
megamaser disks in NGC 1068, NGC 2273, NGC 
4388, and Circinus. This is comparable to the fac- 
tor of 4 accuracy expected for the virial estimation 
method. 

6. The accurate BH masses in the seven megamaser 
galaxies contribute to the observational basis of the 
M — a* relation at the low-mass end. The devia- 
tion from the mean relation of the several accurate 
maser BH masses suggests that the relation may 
not be a single, low-scatter power law as originally 
proposed, which has interesting implications for the 
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APPENDIX 

APPLYING THE VIRIAL ESTIMATION METHOD TO MEGAMASEFl GALAXIES 

The virial estimation method for measuring BH masses in AGNs (e.g. [G reene iDS (pOOl . IKim et al.l d200l . 
Vestergaard & Osmer (2009)) uses the broad-line region (BLR) gas as a dynamical tracer. It is usually only applied 
to Type 1 AGNs, where the BLRs can be directly observed. In this method, one estimates the BH mass with the 
following equation: 

M. = fR ^<° , (Al) 

where /, -Rblr, and a\[ ne have been defined in section 5. Since one can directly detect light from the BLRs in Type 1 
AGNs, ciine can be measured from the broad line spectra and one can use the continnum luminosity ALa(5100A) to 
estimate R B lr via the AL A (5100A) - R B lr correlation (e.g. iKaspi et all (|2000l . I2005T) ). 

In a Type 2 AGN, including the megamaser galaxies we study here, our line-of-sight to the BLR is blocked by heavy 
dust extinction, so one cannot directly measure ALa(5100A) and <7n n c- Instead, one can probe the BLRs in megamaser 
galaxies with polarized scattered light and hard X-rays. Among all the megamaser galaxies with measured BH masses, 
polarized scattered light from the BLRs has been detected in NGC 1068, NGC 4388, NGC 2273, and Circinus (see 
Table 6), and X-ray measurements are also available for these four galaxies. 

We can estimate an ne from the relation (June/VpWHM = 2.09 ± 0.45 from I Woo et al.l |2010), where Vfwhm is the 
FWHM width of polarized scattered light in the broad H Q line (NGC 2273, NGC 4388, and Circinus) or H^ line (NGC 
1068). The major concern with these linewidths is that the observed values may not be the same as the linewidth one 
would measure if the BLRs coul d be observed dire ctly. We identify two effects that can induce such a difference from 
the well-studied case NGC 1068 (Mi ller et al.lll99il) . First, the broad H^ lines in the polarized flux spectra can contain 
a contribution from the narrow H^ lines. Without removing the contribution from narrow lines, the BLR linewidth 
may be underestimated by ~20-30% in NGC 1068. Second, the polarized emission from AGNs may originate from light 
being scattered by electrons with a temperature a few times 10 5 K, which results in significant thermal broadening 
(~50%) of the spectral lines. Because of these two effects, there will be systematic errors in <jii ne if one directly uses the 
observed linewidth of the polarized lines to estimate ciine- Among the four megamaser galaxies we consider here, only 
NGC 1068 has been studied in enough detail to remove these two effects. Luckily, the two effects change the linewidth 
in opposite directions, and hence could offset each other to a certain extent. As in NGC 1068, if no correction is made 
for these two effects, the systematic error will be only ~10%, just slightly larger than the measurement error. Without 
knowing the actual contributions of these two effects for the other three galaxies, we assume that the two effects cancel 
each other to the same extent as in NGC 1068 and use the observed linewidths as the approximations for the intrinsic 
widths. The reader should be aware of this caveat when interpreting the comparisons made. 

In this work, -Rblr was estimated from the "intrinsic" £(2-iokev) °f the nuclear region via the ^(2-iokcv) ~~ -Rblr 
correlation ()Kaspi et al.l I2005T ) . Since three (NGC 1068, NGC 2273, and Circinus) of the four megamaser galaxies 
considered here are Compton-thick (i.e. the X-ray absorbing column density is > 10 24 cm -2 ), we paid particular 
attention to how the £(2-iokcv) were measured. The Compton thick nature of these AGNs is a problem because the 
intrinsic radiation is mostly suppressed and the X-ray spectrum is dominated by the reflected or scattered components. 
It is difficult to measure the actual abso rbing column density and it is very likely that the intrinsic hard X-ray 
luminosity is severely underestimated, e.g. iLevenson et al.l (2006); Bassani etatl (fl999h . Therefore, we excluded those 
measurements that did not consider the Compton-thick nature of these sources and failed to give the absorbing column 
density in the expected range. We mainly considered those measurements from either data with appropriate modeling 
or from observations with instruments capable of directly measuring the transmission components of X-ray above 10 
kev. We took at least two different me asurements f or eac h galaxy from the literature and used the average value to 
calculate Hblr from the correlation in lKaspi et al.l ([2005D . 

Given Rblr and an ne , we estimated the BH masses using Equation Al with the empirically determined (/) = 5.2^'g 
from lWoo et al.1 ([2010). The resultant BH masses were compared with the maser BH masses in section 5. 

APPENDIX 
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TABLE 1 
The Megamaser Sample 





R.A. 


Decl. 


<5RA 


(5DEC 


Vsys 


<5V s ys 


Spectral 


Hubble 


Name 


(J2000) 


(J2000) 


(mas) 


(mas) 


(km s -1 ) 


(km s -1 ) 


Type 


Type 


NGC 1194 


03:03:49.10864 a 


-01:06:13. 4743 a 


0.2 


0.4 


4051 


15 


Sy 1.9 


SAO+ 


NGC 2273 


06:50:08.65620° 


60:50:44.8979° 


10 h 


10 h 


1832 


15 


Sy 2 


SB(r)a 


UGC 3789 


07:19:30.9490 c 


59:21:18.3150 c 


10 h 


10 h 


3262 


15 


Sy 2 


(R)SA(r)ab 


NGC 2960 


09:40:36. 38370 d 


03:34:37.2915 d 


10 h 


10 h 


4945 


15 


Liner 


Sa? 


NGC 4388 


12:25:46. 77914 s 


12:39:43.7516 e 


0.1 


0.3 


2527 


1 


Sy 2 


SA(s)b 


NGC 6264 


16:57:16.12780 f 


27:50:58.5774 f 


0.3 


0.5 


10213 


15 


Sy 2 


S ? 


NGC 6323 


17:13:18.03991 g 


43:46:56.7465 g 


0.2 


0.4 


7848 


10 


Sy 2 


Sab 



Note. — (l)The systemic (recessional) velocities of the galaxies, V S ys, listed here are based on the "optical" velocity convention (i.e. no 
rclativistic corrections are made) , measured with respect to the Local Standard of Rest (LSR). Except NGC 4388, the systemic velocities 
Vsys are obtained from fitting a Keplerian rotation curve to the observed data as described in section 3. The uncertainties 5V S ys given 
here include both the fitting error and a conservative estimate of the systematic error. The fitting error is typically only about 5 km s 
and the systematic error is from possible deviation of the position of the BH fr om (9 x ,6y) = (0 , 0) (see section 3), which is assumed in our 
rotation curve fitting. For NGC 4388, we adopt the observed HI velocity from lLu et al.l 1)20031 ). which has a small but perhaps unrealistic 
error. (2) The positions of UGC 3789, NGC 1194, NGC 6323, and NGC 4388 refer to the location of maser emission determined from our 
VLBI phase-referencing observations. The positions of NGC 2273 and NGC 2960 are determined from K-band VLA A-array observations 
of continuum emission from program AB1230 and maser emission from AB1090, respectively. The maser position for NGC 6264 is derived 
from a phase-referencing observation in the VLBA archival data (project BK114A). (3) The Seyfert types and morphological classifications 
are from NASA/IPAC Extragalactic Database (NED). 

a The position of the maser spot at V op = 4684 kms -1 , where V op is the "optical" velocity of the maser spot relative to the LSR. 
b The position of the radio continuum emission observed at 21867.7 MHz and 21898.9 MHz. 
c The position of the maser spot at V op = 2689 kms -1 . 
d The position of the maser spot at Vop = 4476 kms -1 . 
e The position of the maser spot at V op = 2892 kms -1 . 

1 The average position of the masers with velocities from Vop = 10180 kms -1 to Vop = 10214 kms -1 
g The position of the maser spot at V op = 7861 kms -1 . 

h For all positions derived from VLA observations, we use 10 mas as the actual position error, rather than use the fitted error from the 
VLA data, which are only a few mas for these galaxies. The reason is that the systematic error caused by the imperfect tropospheric model 
of the VLA correlator can be as large as a few to 10s mas. A 10 mas position error usually leads to a 3 — 5 % error in the BH mass for the 
megamasers presented here. Note that although the position for UGC 3789 is derived from a VLBI phase-referencing observation with a 
phase calibrator 2.1° away, the position of this calibrator is derived from a VLA observation. So, we also use 10 mas as the actual position 
accuracy for UGC 3789. 
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TABLE 2 

Observing Parameters 



Experiment 










Synthesized Beam 


Sensitivity 


Observing 


Code 


Date 


Galaxy 


Antcnnas a 


(mas x mas,deg) b 


(mjy) 


Mode c 


BB261 


onntr m c\c 

2UU5 Mar Do 


Ubt 3789 


VLB A, 


GB, EB 


n cr er w n ere on n 

0.55x0.55, 39.0 


1.0 


Self-cal. 


r>r5242r> 


2UU7 JNov 13 


rNGC 1194 


VLBA, 


GB 


2.45x0.36,-7.5 


~ 1.5 


Self-cal. 


T "J T "J O /I O T "\ 


2008 Jan 21 


NGC 1194 


VLB A, 


GB 


1.83x0.44, —10.5 


2.1 


Phase-ref. 


BBzolB 


2UU9 reb 28 


■\T(^~1/~1 00""70 

JNCxL- 2273 


VLBA, 


GB 


n oc\^,r\ on oo a 

U.o9xU.39, —22.4 


0.5 


Self-cal. 


BB2611 1 


2009 Apr 06 


NGC 6264 


VLBA, 


GB, EB 


0.93x0.38,-24.8 


0.5 


Self-cal. 


BB261H 


2009 Apr 18 


NGC 6264 


VLBA, 


GB, EB 


1 no v , n on ~t a t 

1.02x0.29,-14.7 


0.3 


Self-cal. 


T "J T "J O O 1 T~' 

BB231E 


2007 Apr 07 


NGC 6323 


VLBA, 


GB, EB 


nr"Ov,no/i i o c* 

0.58x0.24,-13.6 


0.8 


Phase-ref. 


riri231r 


oha7 A no 

2007 Apr 08 


T\T/~I 1 /"'OOO 

1NGG 6323 


VLBA, 


GB, EB 


0.58x0.27,-14.8 


0.8 


Phase-ref. 




nnnT A nn 

2007 Apr 09 


NGG 6323 


VLBA 


GB EB 


n no w n oo tr o 

0.92x0.32,-5.8 


0.8 


Phase-ref. 


BB231H 


2007 Apr 29 


NGC 6323 


vlba| 


GB| EB 


/~\ o n v > n o o 10 ft 

0.90x0.28,-13.5 


1.3 


Phase-ref. 


jtSrs242.b 


2UU8 Apr 13 


■\T r~\C* go no 

JNCrC o323 


VLBA, 


GB 


U.95XU.44,— 9.4 


0.8 


T>T „„ r 

jrnasc-rci. 


RR949F, 


900K Anr U 


IN V_xv_' UOZO 


VLBA, 


GB 


QOvO 4Q — Q 9 


0.8 


1 * n *"i on rot 


BB242G 


2008 Apr 16 


NGC 6323 


VLBA, 


GB 


0.79x0.45,-10.2 


1.0 


Phase-ref. 


BB242H 


2008 Apr 19 


NGC 6323 


VLBA, 


GB 


0.80x0.45,-4.2 


0.3 


Self-cal. 


BB242J 


2008 May 23 


NGC 6323 


VLBA, 


GB, EB 


0.58x0.21,-8.1 


1.7 


Phase-ref. 


BB242M 


2009 Jan 11 


NGC 6323 


VLBA, 


GB, EB 


0.52x0.29,-17.8 


0.4 


Self-cal. 


BB242R 


2009 Apr 17 


NGC 6323 


VLBA, 


GB, EB 


0.87x0.45,12.2 


0.5 


Self-cal. 


BB242S 


2009 Apr 19 


NGC 6323 


VLBA, 


GB, EB 


0.57x0.26,-21.2 


0.4 


Self-cal. 


BB242T 


2009 Apr 25 


NGC 6323 


VLBA, 


GB, EB 


0.75x0.23,-19.0 


0.8 


Self-cal. 


BB248 


2009 Mar 07 


NGC 2960 


VLBA, 


GB, EB 


0.97x0.45, 5.3 


2.5 


Self-cal. 


BB261C 


2009 Mar 20 


NGC 2960 


VLBA, 


GB, EB 


1.11x0.48,2.6 


1.1 


Self-cal. 


BB261D 


2009 Mar 23 


NGC 2960 


VLBA, 


GB, EB 


1.26x0.44, -2.1 


1.3 


Self-cal. 


BB184C 


2006 Mar 26 


NGC 4388 


VLBA, 


GB, EB 


1.23x0.32, -9.27 


1.8 


Phase-ref. 



a VLBA: Very Long Baseline Array; GB: The Green Bank Telescope of NRAO; EB: Max-Planck-Institut fur Radioastronomie 100 m 
antenna in Effelsberg, Germany. 

b Except for program BB184C, this column shows the average FWHM beam size and position angle (PA; measured east of north) at the 
frequency of systemic masers. For BB184C, the FWHM and PA is measured at the frequency of red-shifted mascrs because no systemic 
maser is detected (FWHM and PA differ slightly at different frequencies). 

c "Self-cal." means that the observation was conducted in the "self-calibration" mode and "Phase-ref." means that we used the "phase- 
referencing" mode of observation. 



TABLE 3 

Sample data for NGC 6264 



V a 

V op 


RA b 


5RA b 


Decl. b 


<5DEC b 


F c 




(km s- 1 ) 


(mas) 


(mas) 


(mas) 


(mas) 


(mJy/B) 


(mJy/B) 


10918.33 


0.399 


0.008 


-0.021 


0.018 


6.0 


0.3 


10914.71 


0.407 


0.010 


-0.001 


0.021 


4.8 


0.3 


10911.09 


0.401 


0.005 


0.000 


0.011 


10.2 


0.4 


10907.47 


0.397 


0.013 


0.041 


0.029 


3.7 


0.3 


10903.85 


0.395 


0.011 


-0.050 


0.027 


3.8 


0.4 


10900.23 


0.388 


0.009 


0.044 


0.022 


4.6 


0.3 


10885.74 


0.369 


0.019 


0.053 


0.038 


2.6 


0.4 


10871.26 


0.400 


0.018 


0.061 


0.037 


2.9 


0.3 


10856.77 


0.467 


0.011 


-0.038 


0.029 


3.8 


0.3 


10853.15 


0.484 


0.005 


-0.038 


0.011 


9.7 


0.4 


10849.53 


0.490 


0.005 


-0.006 


0.011 


9.4 


0.3 


10845.91 


0.492 


0.002 


-0.016 


0.006 


18.3 


0.4 


10842.29 


0.494 


0.002 


-0.005 


0.005 


24.4 


0.4 


10838.67 


0.505 


0.004 


-0.010 


0.009 


11.7 


0.4 


10835.05 


0.543 


0.012 


-0.094 


0.032 


3.5 


0.4 


10831.43 


0.542 


0.017 


-0.039 


0.031 


3.1 


0.3 


10827.81 


0.518 


0.012 


-0.101 


0.027 


3.9 


0.3 


10824.19 


0.533 


0.007 


-0.040 


0.016 


6.9 


0.4 


10820.56 


0.493 


0.012 


-0.041 


0.026 


3.8 


0.3 


Note. - 


Sample of data for NGC 6264. The entirety of data for all galaxies 


is available in 


the electronic version. 





a Velocity referenced to the LSR and using the optical definition (no relativistic corrections). 

b East-west and north-south position offsets and uncertainties measured relative to the average position of the systemic masers in the 
VLBI map (Figure 2). Position uncertainties reflect fitted random errors only. In NGC 1194 and NGC 2960, there may be additional 
uncertainties caused by the poorer tropospheric delay calibrations due to their low declinations (low elevations during the observations). 
c Fitted peak intensity and its uncertainty in mjy beam -1 . 
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TABLE 4 

The BH Masses and Basic Properties of the Maser Disks 





Dist. 


BH mass 


Disk Size 


P.A. 


Incl. 


r c 


P0 




T a go 


Rinf 


Name 


(Mpc) 


(10 7 M ) 


(PC) 


(°) 


(°) 


(pc) 


(M Q pc" 3 ) 




(yr) 


(arcsec) 


NGC 1194 


52 


6.4±0.3 


0.53-1.30 


157 


85 


0.250 


1.2xl0 9 


12 


>1.0xl0 1(1 


0.033 


NGC 2273 


26 


0.76±0.04 


0.028-0.085 


153 


84 


0.015 


6.0X10 11 


0.05 


<2.2xl0 6 


0.010 


UGC 3789 


50 


1.12±0.05 


0.09-0.32 


41 


> 88 


0.024 


2.0X10 11 


0.12 


<2.7xl0 7 


0.010 


NGC 2960 


71 


1.14±0.05 


0.13-0.36 


-131 


89 


0.055 


1.8xl0 10 


0.5 


<9.3xl0 8 


0.005 


NGC 4388 


19 


0.84±0.02 


0.24-0.29 


107 




0.090 


3.3xl0 9 


0.9 


<6.6xl0 9 


0.034 


NGC 6264 


136 


2.84±0.04 


0.23-0.78 


-85 


90 


0.083 


1.3xl0 10 


1.3 


<8.0xl0 9 


0.012 


NGC 6323 


105 


0.93±0.01 


0.13-0.30 


10 


89 


0.046 


2.3xl0 10 


0.3 


<4.7xl0 8 


0.003 



Note. — Col(l): Galaxy name; Col(2): Distances we adopt from NED; Col(3): Black hole masses measured in our project. The mass 
uncertainty here only includes errors caused by source position uncertainty and from fitting a Keplerian rotation curve for a given distance. 
The actual uncertainty is limited by the error of the Hubble Constant Ho (11%); Col(4): Sizes of the inner and outer edge of the maser 
disks; Col(5): Position angle (PA) of the disk plane measured east of north. PA equals zero when the blueshifted side of the disk plane has 
zero East offset and positive North offset; Col(6): Inclination of the maser disk. Note that the inclination of NGC 4388 is unconstrained 
because we did not detect systemic masers; Col(7): The core radius of the Plummer cluster in parsecs. For NGC 1194, NGC 2273, NGC 
2960, and NGC 4388, the core radii are derived from Equation 12. For UGC 3789, NGC 6264, and NGC 6323, the radii are derived from 
the Plummer rotation curve fitting; Col(8): The central mass density of the Plummer cluster po = 3Moo/47rrj?. Here, Moo is obtained 
from the Plummer rotation curve fitting ; Col(9): The maximum stellar mass of the Plummer cluster below which the cluster will not 
evaporate in less than 10 Gyr. In all cases except NGC 1194, a cluster of neutron stars can be directly ruled out because m max is less 
than Ri 1.4 Mq; Col(10): Lifetime (T age ) of a cluster. The values shown here are limited by the collision timescale, which is the maximum 
lifetime of the cluster composed of either main-sequence stars, brown dwarfs, white dwarfs, or neutro ns stars; Col ( ll) T he radius of the 
gravitational sphere of influence for the maser BHs in arcsec. We calculate Ri„f using Equation 1 of [Barth ct al. (2003) with the bulge 
velocity dispersion measurements from Greene et al. (20l3). 



TABLE 5 

Upper limit on Continnum Emission from Megamaser Galaxies 



Galaxy 


Vcenter 


VRango 


ha 


Project Code 




(km s^ 1 ) 


(km S" 1 ) 


(mjy) 




NGC 1194 


3568 


3489 - 3660 


< 0.34 


BB242B 


NGC 2273 


2800 


2868 - 2909 


< 0.18 


BB261B 


NGC 2960 


5315 


5296 - 5354 


< 0.14 


BB261C, BB261D 


NGC 4388 


2600 


2881 - 2966 


< 0.42 


BB184C 


NGC 6264 


9710 


9726 - 9776 


< 0.16 


BB261F, BB261H 


NGC 6323 


8100, 7650 


8005 - 8200, 7556 - 7750 


< 0.08 


BB242M, BB242R, BB242S, BB242T 



Note. — Col(l): Name of the galaxy; Col(2):The central Optical-LSR velocity of the bands used to search for continuum emission; 
Col(3) The velocity range corresponding to the selected channels. Note that the channels selected are chosen to be free of any maser lines 
except for NGC 2960. For this galaxy, we don't have line-free channels, and we averaged the bands centered at the systemic velocity of 
the galaxy and searched for continuum emission offset from the systemic maser emission; Col(4) The 2a detection limit of the continuum 
emission; Col(5) The data used for continuum detection. 



TABLE 6 

Comparison of Maser BH Mass with Mass from Virial Estimation 



Galaxy 


maser BH 


virial BH 




Vfwhm 


Reference 


L x (2-10 kev) 


Reference 




(10 6 Mq) 


(10 6 M ) 


(light-days) 


(km s" 1 ) 




(10 42 ergs s" 1 ) 




NGC 1068 


8.6±0.3 


9.0±6.6 


4.6±3.0 


2900 


1 


6.5, 2.6 


5, 6 


NGC 4388 


8.5±0.2 


7.1±4.9 


2.0±1.2 


3900 


2 


0.8, 0.9,1.9,1.0 


7,8,9,10 


NGC 2273 


7.6±0.1 


4.3±2.8 


2.2±1.2 


2900 


3 


1.0, 1.7 


11,12 


Circinus 


1.7±0.3 


4.8±3.2 


1.9±1.1 


3300 


4 


1.1, 1.0, 1.2 


13,14,15 



Note. — Col( l ): Ga laxy name; Col(2) BH mass measur e d from the megamaser technique. The BH mass of NGC 1068 is taken from 
ILodato fc Bertinl (I2003T) and Circinus from IGreenhill et al.l J2003T) . The BH masses of NGC 4388 and NGC 2273 are from this paper. 
Col(3): BH mass measured fro m the virial estim ation method. Col(4): The size of the Broad Line Region (BLR) calculated using the 
L(2— lOkev) ~ R-SLii correlation ( IKaspi et al.li2005l ). Col(5): The full width at half maximum (FWHM) of the observed broad line (Hg for 
NGC 1068 and H a for NGC 4388, NGC 2273, and Circinus). In NGC 2273, only the full width at zero intensity (FWZI) is given, so we 
estimate the FWHM "by eye". We assume the measurement uncertainty is 200 km s — 1 in all cases. Col(6) Reference for the linewidth 
measurement. Col(7) Intrinsic hard X-ray (2-10 kev) luminosity. We used the average luminosity to calculate the size of BLR. The error 
of the luminosity is taken to be the standard deviation of the adopted luminosities. Note that the distances used to calculate the intrinsic 
luminosities are 15.4, 19.0, 26.0, and 4.0 Mpc for NGC 1068, NGC 2273, NGC 4388, and Circinus, respectively. Col. (8) Reference for the 

adopted X-ray lu minosities. 

Refere nce s - 1. IMiller et al l |[T99"iT) 2. IHo et al.l ([19971) 3. rivToran et al.l (127)001) 4. IQlival fTSgSl) 5. ILevenson et al I 12006T) 6. IQgle et al. 
(2003) 7. ICappi et alj (12009) 8. lAkvlas fc*G corgantopoulos ( 2 0091) 9. IBassani et aT. (1999) 10. IForster et al.l 119991) 11. Tcrashim a et al, 
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Fig. 1. — Characteristic H2O maser spectra. The x-axis shows LSR velocities based on optical definition. Flux densities of masers can 
vary significantly, so the spectra shown here are just representative for particular epochs: January 13 2008 for NGC 1194; February 21 
2009 for NGC 2273; April 2 2009 for NGC 2960 (Mrk 1419); November 30 2005 for NGC 4388; March 31 2009 for NGC 6264; and April 6 
2000 for NGC 6323. 
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Fig. 2. — VLBI maps for the seven 22 GHz H2O masers megamasers analyzed. The maps are color-coded to indicate redshifted, blueshifted, 
and systemic masers, where the "systemic" masers refer to the maser components having recessional velocities close to the systemic velocity 
of the galaxy. Except NGC 4388, maser distributions are plotted relative to the average position of the systemic masers. For NGC 4388, 
in which the systemic masers are not detected, we plot the maser distribution relative to the dynamical center determined by fitting the 
high velocity features with a Keplcrian rotation curve. 
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Fig. 3.— Maser distributions (top panels) and rotation curves (bottom panels) for NGC 1194 , NGC 2273, UGC 3789, and NGC 2960. 
The maser distribution has been rotated to horizontal to show the scatter in the maser positions and the offset of the systemic masers from 
the plane defined by high-velocity masers more clearly. The coordinate system is chosen to place the centroid of the high-velocity maser 
disk (blue and red points) at 8 y = and the centroid of the systemic masers (green points) at X = 0. The axes for the maps show relative 
position in milliarcseconds, and North (N) and east (E) are indicated by directional arrows on each map. The bottom panel for each galaxy 
shows the rotation curves of the redshifted and blueshifted masers (red and blue points on the curves) plotted with the best-fit Keplerian 
(solid curve) and Plummer (dotted curve) rotation curves. The velocities shown in the figure are the LSR velocities after the special and 
general relativistic corrections. The residuals (data minus Keplerian curve in red and blue; data minus Plummer curve in black) are in the 
bottom part of each figure. Note that we plot the rotation curve with the impact parameter 9 (mas) as the ordinate and rotation speed \v\ 
(km s~ 1 ) as the abscissa for the convenience of fitting. 
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Fig. 4.— Maser distributions (Top panel) and rotation curves (Bottom panel) for NGC 4388 , NGC 6264, and NGC 6323. Please refer 
to the caption of Figure 3 for the description of this figure. 



